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Abstract: We review the equivalence of maximal center gauge fixing to the problem 
of finding the best fit, to a given lattice gauge field, by a thin vortex configuration. 
This fit is necessarily worst at the location of P-plaquettes. We then compare the fits 
achieved in Gribov copies generated by (i) over-relaxation; (ii) over-relaxation after 
Landau gauge preconditioning; and (iii) simulated annealing. Simulated annealing 
yields the best fit if all links on the lattice are included, but the situation changes if we 
consider only the lattice volume exterior to P-plaquettes. In this exterior region, the 
fit is best for Gribov copies generated by over-relaxation, and worst for Gribov copies 
generated after Landau gauge preconditioning. The two fitting criteria (including or 
not including the P-plaquettes) yield string tensions differing by —34% to +20% 
respectively, relative to the full string tension. Our usual procedure ("quenched 
minimization") seems to be a compromise between these criteria, and yields string 
tensions at an intermediate value close to the full string tension. 
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1. Introduction 

Direct maximal center gauge is the common name for lattice Landau gauge in the 
adjoint representation. For the SU(2) group, adjoint links Ua in this gauge satisfy 
the condition 

TTY,L,{UA^{x)-UAf.{x-fl)) =0 (1.1) 

at every point (the {Li} are the SU(2) group generators in the adjoint representation) 
leaving a residual Z2 symmetry. This gauge is used in studies of confinement, in 
particular to locate center vortices in thermalized lattice configurations. 

As in ordinary lattice Landau gauge, there are many points on the gauge orbit, 
not related by any residual symmetry, which satisfy the local gauge condition (|1.1| ). 
These configurations are known as Gribov copies. Difficulties associated with Gribov 
copies in maximal center gauge have been noted in recent work by Bornyakov et 
al. [1,2], and previously by Kovacs and Tomboulis [3]. In this article we discuss 
some issues that we believe are relevant to the problem of generating a set of Gribov 
copies in maximal center gauge, and to choosing the "best" Gribov copy among 
the set. Before proceeding, we should note that the Gribov problem can be avoided 
altogether by using the Laplacian center gauge [4] . On the other hand, the Laplacian 
version does have certain unattractive features, notably the lack of scaling of the P- 
vortex density [5] and the lack of "precocious linearity" [6]. While these features 
of Laplacian center gauge are by no means fatal, we think it worthwhile to further 
explore the Gribov copy issue in the older proposal of maximal center gauge. 
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2. Gauge Fixing as a "Best Fit" Procedure 

Gauge fixing always has a flavor of arbitrariness, and claims for a privileged gauge 
are typically regarded with suspicion. To understand the rationale for maximal 
center gauge, it is essential to understand its function as a flt of a given lattice 
configuration by a singular pure gauge (thin vortex) field. This insight, which we 
will now elaborate, is due to Engelhardt and Reinhardt in ref. [7]. 

Imagine running a Monte Carlo simulation of lattice SU(2) gauge theory at high 
/3, and printing out the values of link variables U^{x) in some thermalized configura- 
tion. At a glance, these values would look like random numbers, but of course this 
impression is deceptive. Locally the link variables are only small fluctuations around 
a pure gauge configuration. Suppose we then ask for the pure gauge configuration 
g{x)g^{x + ft) which is closest, in lattice configuration space, to the given lattice 
gauge field U^j,{x). Using the standard metric on the SU(2) group manifold, this is 
equivalent to asking for the gauge transformation g{x) which minimizes the square 
distance in configuration space 

= Tr - g{x)g\x + mU,{x) - g{x)g\x + fi))^' 

= 2 E(l^[^2] - Tr[g\x)U,{x)g{x + fi)]) (2.1) 

X,IJ. 

This quantity is minimized by finding the gauge-transformed configuration 

^U^{x) = g\x)U^{x)g{x + fi) (2.2) 

which maximizes 

Rian = Y.T^r[m^{x)] (2.3) 

Xy/J, 

From this we conclude that the problem of finding the pure-gauge configuration 
closest to a given lattice gauge field is completely equivalent to the problem of fixing 
that lattice gauge field to the Landau gauge. 

We next generahze this idea shghtly, and allow for Z2 dislocations in the gauge 
transformation. This means fitting the lattice configuration by a slightly more general 
form 

U;'''-{x)^g{x)Z,{x)g\x + fi) (2.4) 

where Z^{x) ~ ±1. This is a thin center vortex configuration, generated by a sin- 
gular gauge transformation. Note that C/^°^ becomes a continuous pure gauge in the 
adjoint representation, which is blind to the Zn{x) factor. This motivates a two-step 
fitting procedure: First, we determine g[x) up to a residual Z2 transformation, by 
minimizing the square distance d\ in configuration space between t/^(a;) and U^°'^{x) 
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in the adjoint representation: 
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^Tr[(f/A^(x) - gA{x)g\{x + f^)){UA^J,{x) - gA{x)gli{x + p,)^ 

X,fl 

2j2(Ti[h]-Ti[g\{x)UA^ix)gA{x + il)] 

X,fl 

2 5] (4- TT[g\x)U,{x)g{x + fi)] 



(2.5) 



x,ii 



The distance function in this case is determined by the metric on the SU {2)/Z2 group 
manifold. Minimizing d\ is equivalent to fixing to direct maximal center gauge, which 
seeks the largest value of 



Rdn.c = T.\'^T[m,{x)] 



(2.6) 



X,IM 



The condition that R is stationary leads to the adjoint Landau gauge condition (|1.1| ). 
Having determined g{x), we then find Z^{x) at each link by minimizing 

llix) = Tt[{U,{x) - gix)Z^ix)g^ix + mU,{x) - gix)Z^ix)g^ (x + /i))t] (2.7) 

which is easily seen to require 



Z^,{x) 



signTrPf/^(a;)] 



(2. 



This is the center projection prescription. 

In this way we have shown that maximal center gauge, together with the rule 
( P^.8D for center projection, is equivalent to finding a "best fit" of a given lattice 
configuration U^{x) by a thin vortex configuration f/^°'"(x). While this (rather trivial) 
derivation may be novel, the essential point has been made previously in ref. [7]. 



2.1 Center Dominance 

The gauge-transformed lattice configuration in maximal center gauge can be written 
m^{x) = Z^(x)e'^'^(") (Tre'^^(-) > 0) (2.9) 
with the original lattice configuration 

U,{x) = gix)Z^{x)e'^-^^^g\x + fi) (2.10) 

where g{x), Z^{x) are determined by the procedure just described. Our claim, elab- 
orated in refs. [6,8], is that the confining properties of U^{x) are entirely encoded in 
the Z2 link variables Z^{x)\ the variables A^{x), which are typically of order 1/(3, 
are responsible for short range effects such as the Coulomb force law. Two necessary, 
but not sufficient, conditions for this claim to hold true follows: 
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1. Set Afj^{x) = in ( p.lO|) , which sends f/^(x) — > f/™'"(x). The asymptotic string 



tension in the t/™'" configuration should match the asymptotic tension of the 
full configuration, i.e. "center-projected" Wilson loops 

Z{C) = ^Tr[?7"°"(C)] (2.11) 

have the same string tension as full Wilson loops (where Z{C) denotes a product 
of center-projected links around loop C). This property is known as "Center 
Dominance." 



2. Dropping the Z^{x) variables in (|2.10|) , which can be interpreted as remov- 



ing center vortices from the system, the Wilson loops constructed from the 
remaining degrees of freedom 

W{C) = Tr[ n e'^^] 

linksSiC 

= Z{C)T:t[U{C)] (2.12) 

should have zero string tension. This test was carried out by de Forcrand and 
D'Elia in ref. [9]. 

The reason that these properties are not sufficient for our purposes is that confin- 
ing fiuctuations might be hidden in the Z^{x) variables along with a lot of short-range 
physics. If that is the case, then we might not gain much from the truncation of 
to Z^. For example, even a naive projection Z^{x) = signTr[?7^(x)] of the original 
configuration (with no gauge fixing) has the center dominance property, as is easily 
verified by simple group-theoretic arguments [6,10]. But this naive projection also 
reproduces — exactly! — the Coulomb potential at short distance scales, which is 
certainly due to gaussian field fiuctuations rather than a vortex mechanism. Our 
objective is to strip away such short-range effects, and to isolate as far as possible 
the degrees of freedom which are solely responsible for infrared physics. 

2.2 Where Fitting Fails 

The square deviation at each link, away from a pure gauge in the adjoint represen- 
tation, is proportional to 

^(^) = ^Tr[(t/A;.(x) ~ gA{x)g\{x + fi)){UA^,{x) - gA{x)g\{x + fi)^ 

= I - \{Tr[m,ix)]y (2.13) 

This quantity is generally small, of order 1/(3, at a local minimum of d\. The 
exception is for some links (at least one) belonging to each P-plaquette. For these 
links we must have a very poor fit, with 5^(x) ~ 0(1) regardless of (3. We recall 
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that a plaquette p is a P-plaquette iff Z{p) = —1, and tliat P-plaquettes belong to 
P-vortices. 

The fit to a thin vortex has to be bad at P-plaquettes for the following reason: 
At large P, we generally have (even at P-plaquettes) 

Itt[U{p)] = 1 - O(^) (2.14) 

On the other hand 

Tr[U{p)] = Z{p)Tr[ J] e'^^^"^] (2-15) 

linksSip 

But if Z{p) = —1, equations ( |2.14| ) and ( |2.15| ) imply that Af^{x) is 0(1) rather than 
0(1/ (3), on one or more links in the plaquette p. The magnitude of A^{x) is a measure 
of the goodness-of-fit, and in fact S"^ depends only on this variable. A perfect fit on 
a link corresponds to = 0, while ~ 0(1) is a very bad fit. 

In ordinary Landau gauge, S'^{x) can in principle be small on every link. In 
contrast, in maximal center gauge, we see that this property can only hold if there 
are no P-vortices found on the lattice. 



3. Quenched vs. Annealed Gribov Copies 

There is no known method for finding the global minimum of the distance function 
d\, but two methods have been employed to generate local minima, i.e. Gribov copies, 
satisfying the local condition (|1.1|). These are the methods of simulated quenching, 
and simulated annealing. The terms derive from the fact that minimizing d\ is 
equivalent to finding the zero-temperature ground state of an analog "spin-glass" 
system 

H = -J2^'^[9A{x)UAt,{x)gA{x + fi) 

X,fJ, 

= (TrM(x)f/^(x)5f(X + /i)j) +const (3.1) 



where gA^x) is the dynamical S0(3) group-valued "spin" variable, and Ua is a set of 
fixed, stochastic, nearest-neighbor couplings. 

An obvious approach to finding the zero-temperature ground state is to cool 
the system, either gradually ("annealing") or suddenly ("quenching"). Quenching 
corresponds to placing the system in contact with a reservoir at zero temperature, 
and implies that only changes in g{x) which lower the energy density can be accepted. 
In practice, quenching is implemented by the over-relaxation method, as described in 
ref . [8] . Simulated anneahng is achieved by applying Metropolis updates to the analog 
spin system, and gradually lowering the temperature variable from some initial value 
to zero [1]. Neither of these methods obtains the true minimum of H; each generates 
a set of Gribov copies in maximal center gauge. 
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A variation of the quenching approach was studied by Kovacs and Tombouhs 
[3]. Instead of applying over-relaxation to a random point on the gauge orbit of 
a thermalized configuration, they first fixed the configuration to ordinary lattice 
Landau gauge, and then apphed over-relaxation to fix to maximal center gauge. We 
will refer to the quenching procedure via over-relaxation described in [8] as "OR", 
the Kovacs- Tomboulis variation of this procedure as "KT" , and simulated annealing, 
studied in maximal center gauge by Bornyakov et al. [1], as "SA". 

In Table |l| we display the average deviation 5^ per link 

5' = 4^E^» (3.2) 

for Gribov copies generated by the KT, SA, and OR methods {4V is the number 
of links on the lattice). These methods were applied to thermalized lattice SU(2) 
configurations generated at the f3 values and lattice sizes shown below. Only one 
Gribov copy was generated, in each method, for each thermalized lattice; there was 
no attempt to select the "best" copy out of a set. We also display, in Fig. 0, the 
fractional deviation 

6^-6\OR) 



f 



S^{OR) 



(3.3) 



of S"^ away from the OR result obtained in each method. We see that, according to 
the criterion of minimizing the average S"^, simulated annealing gives the best fit to 
a thin vortex, the KT method is a little worse, and OR is the worst of the three, as 
already noted in ref. [1]. 





lattice size 


Kovacs- Tomboulis 


sim anneal 


pure over-relax 


2.1 


10^ 


0.2797(3) 


0.2736(3) 


0.2805(4) 


2.2 


12^ 


0.2631(2) 


0.2591(2) 


0.2657(2) 


2.3 


16^ 


0.2440(1) 


0.2406(1) 


0.2467(1) 


2.4 


20^ 


0.2232(2) 


0.2206(2) 


0.2260(2) 



Table 1: Values for the average 6'^, for Gribov copies generated by three different methods. 

There is something a little peculiar about this ordering, however. The OR Gribov 
copies are known to have reasonably good center dominance properties [8,14]. Gribov 
copies generated by the KT method result in a vanishing projected string tension [3], 
while the SA procedure yields projected tensions at roughly 66% of the full string 
tension [1]. One might have expected that the projected string tension would be 
correlated directly with 5^. Instead, what seems to happen is that the projected 
string tension first drops to zero, as 6"^ falls below the OR result in KT copies, and 
then rises again to 66% of the full string tension as 6"^ falls further in SA copies. 

This raises the question of whether there is some other aspect of the fit to a 
thin vortex which is better correlated with the projected string tension. We have 
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Figure 1: Percent deviation (/ x 100) of 5^ from the over-relaxation value, for the simulated 
annealing and Kovacs-Tomboulis methods. Simulated annealing is best at minimizing 6'^. 

already noted that the fit is necessarily very bad at the location of P-plaquettes. This 
suggests just ignoring the quality of fit at the P-plaquettes, where it is guaranteed to 
be bad, and concentrating on the fit in the lattice volume exterior to P-plaquettes. 
We therefore calculate the average square deviation in the exterior region 

Slt = ^T.^» (3-4) 

^^ext ext 

where N^xt is the number of all links not belonging to P-vortices, and the sum runs 
over the set of these external links. The quantity is a measure of the quality of 
fit to a thin vortex in this exterior region. 

Table | displays the values of 6"^^^, and Fig. ||the corresponding fractional devi- 
ation 

_ 51, - 5U0R) 

- 6U0R) ^'-'^ 

in the external region. Here the order is quite different than in Table ^ and Fig. |l[ 
This time, the best fit is achieved by Gribov copies generated by quenching (OR), 
and these copies also have the largest projected string tension. The SA fit is a little 
worse, with the projected string tension ^ 34% lower than the full string tension, 
while KT copies have by far the worst fit, and also have negligible projected string 
tension. From this data, it appears that the projected string tension correlates much 
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better with the fit in the external region, based on 5^^^, than with the overall fit 5^ 
which includes the P- vortex volume. 





lattice size 


Kovacs-Tomboulis 


sim anneal 


pure over-relax 


2.1 


10^ 


0.2062(4) 


0.1951(3) 


0.1936(3) 


2.2 


12^ 


0.2081(4) 


0.1970(3) 


0.1944(3) 


2.3 


16^ 


0.2086(1) 


0.1975(2) 


0.1948(1) 


2.4 


20^ 


0.2043(1) 


0.1950(1) 


0.1927(1) 



Table 2: Values for ^g^,^. 
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Figure 2: Same as Fig. ||, but for the exterior region. Pure over-relaxation is best at 
minimizing (5g^.j. 

The results just obtained suggest selecting among Gribov copies generated by 
quenching, on the basis of the smallest average 51^^. But this procedure is also not a 
great success as regards center dominance. We have found that choosing the lowest 
(5g2.j out of a set of quenched copies leads to projected string tensions which are ~ 20% 
higher {(3 = 2.4) than the full string tension, as seen, e.g., from the projected Creutz 
ratio Xcp(5, 5) shown in Fig. ^. 



4. Quenched Minimization 

If the criterion for the best gauge copy is the best average fit over all links, including 
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|3=2.4 



Bali et al. 

20'^ lattice ^ 



5 10 15 20 

Ncopy 

Figure 3: Xcp(5,5) vs. the number Ncopy of Gribov copies generated per configuration, at 
P = 2.4. Tlie Creutz ratio is evaluated in the copy with the best exterior fit, minimizing 
(5g^^. Solid and dotted lines are the asymptotic string tension and errorbars, respectively, 
on the unprojected lattice [11]. 

P-vortices, then the projected string tension is about 34% lower than the full string 
tension, as shown by Bornyakov et al. in ref. [1]. On the other hand, if the criterion 
for best fit is the fit to a thin vortex in the region exterior to P-vortices, then the 
projected string tension comes out around 20% too high (/3 = 2.4). A case can be 
made for either fitting criterion, but the method which we have used in the past now 
appears to be something of a compromise between the two. We begin by quenching 
the analog spin system ( |3.1| ), starting from each of a set of random points on the 
gauge orbit, to generate a set of Gribov copies satisfying (|1.1|). This procedure, as 
seen in the previous section, tends to emphasize the exterior fit 6"^^^ to a thin vortex. 
Among the gauge copies generated by quenching, we then select the copy with the 
best overall fit; i.e. the minimum value of 5^. 

In the absence of a compelling reason to prefer the 5^ or S^^^ criterion, the 
"quenched minimum" prescription seems a priori as good as any. It is, at least, 
a perfectly well-defined procedure. We have argued in the past, on the basis of a 
number of empirical tests, that center projection applied to quenched minima locate 
physical objects, namely the center vortices. The results of those tests have not 
changed any, and are worth summarizing once again: 

1. Let Wn{C) denote the expectation value of a Wilson loop constructed from 
unprojected links, evaluated in the subensemble of configurations in which n 



un 
un" 
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P-vortices, on the projected lattice, pierce the minimal area of loop C . It is 
found that 

as expected if P-vortices locate center vortices in the unprojected configurations 



2. If a vortex is inserted "by hand" (via a singular gauge transformation) into a 
thermalized lattice, then the set of P-vortices on the projected lattice includes 
the inserted vortex [12]. 

3. Using information about P- vortex location to remove center vortices from the 
unprojected configuration, one finds that the confining and chiral symmetry 
breaking properties of the configuration are also removed, and the topological 
charge goes to zero [9]. 

4. The density of P-vortices, and therefore the density of vortices identified on the 
unprojected lattice, scales correctly according to the renormalization group, as 
first noted in [13] (see also the later results in [6,8]). 

5. Center dominance: The string tension of the projected lattice agrees fairly well 
with the asympotic string tension of the unprojected lattice [8,14]. 

6. Precocious linearity: Projected Creutz ratios Xcp{I, I) vary only a little with 
J, and there is no Coulombic force at small distances. This indicates that the 
projected degrees of freedom are not mixed up with short-range physics [6]. 

Despite these apparent successes, there is still a serious objection that can be 
raised to the quenched minimization approach, in view of the findings of Bornyakov 
et al. [1]. The problem is that any local minimum of d\ can be reached by quenching, 
if the starting configuration is within the "basin of attraction" of that minimum on 
the gauge orbit. Thus, if we sample enough random configurations on the gauge 
orbit by the OR approach, eventually minima obtained by the SA method would 
be reached, and the projected string tension would drop well below the full string 
tension. 

The answer to this objection is based on the fact that the projected string tension 
is found to converge rapidly, with the number of Gribov copies generated, to a value 
in good agreement with the full string tension. Moreover, the convergence improves 
as the lattice size increases [14]. The implication is that the measure of SA copies 
must be negligible compared to the measure of OR copies, at least at large volume. 
Both the volume and copy-number dependence of the projected Creutz ratio Xcp(5, 5) 
at /3 = 2.5 are illustrated in Fig. ^ taken from ref. [14]. Similar examples of other 
Creutz ratios, at various (3 values, can also be found in that reference. The indications 
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are that on an infinite lattice there is convergence to a result well above the SA 
value. The probability of a random configuration evolving, under quenching, to an 
SA minimum seems likely to go to zero in the infinite volume limit (although this 
conjecture deserves further investigation). 
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Figure 4: Xcp(5, 5) vs. the number Ncopy of Gribov copies generated per configuration, 
obtained at various lattice volumes. The quenched minimization procedure is used. 

A last point which we would like to make, in connection with quenched Gribov 
copies, is that the projected lattices of different copies seem to be quite well correlated 
at distances beyond one fermi, which is the width of a center vortex.^ To illustrate 
this fact, we generate two Gribov copies for each thermalized lattice by applying the 
over-relaxation algorithm both to the original thermalized lattice, and to a random 
gauge copy of that lattice. The two gauge-fixed lattices are center projected to 
obtain two projected lattices, denoted Z^[x) and Z^{x). We then calculate projected 
Creutz ratios Xprod{I, I) from the product Wilson loops 

WprodiC) = (Z^(C)Z^(C)) (4.2) 

If the two projected configurations were perfectly correlated, or if the A and B loops 
differed only by perimeter effects, then we would find 

Xprodil, /) = (strong correlation) (4-3) 

^The width of center vortices can be determined in three different ways: from the falloff of 
Wi/Wo [14], from the vortex free energy in a finite volume [15], and from the adjoint string-breaking 
scale [16]. The three determinations are in rough agreement. 
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At the other extreme, if there were no correlation at all between the projected con- 
figurations, we would find 

Xprodil, I) = 2xa(/, I) = 2xi?(/, I) (no correlation) (4.4) 

where x^(/,/) = xb{I,I) are the projected Creutz ratios obtained from Z^(C) or 
Z^{C) loops separately. 

Creutz Ratios at p=2.2 from Products of Z2 Loops 

0.5 I > > > > > > > 1 



0.4 



0.3 - 




-0.1 I ' ' ' ' ' ' ' 1 

1 1.5 2 2.5 3 3.5 4 4.5 5 

I 



Figure 5: Xprod{I, /) at /3 = 2.2. 

Figures || and || show our results for Xprod{I, /) at /5 = 2.2 and P = 2.3, which 
indicate a tendency towards strong correlation of the two projected Gribov copies in 
the infrared. 

5. Conclusions 

We have pointed out that there are at least two reasonable criteria for selecting 
among Gribov copies in maximal center gauge. One can select the copy which is a 
best fit to a thin vortex over the entire lattice volume, including links in P-vortices, 
and this leads to projected string tensions which are some 34% below the full string 
tension. Alternatively, given that the fit is bound to fail at P-plaquettes, one can 
select on the basis of the best fit in the lattice volume exterior to P-plaquettes, and 
in fact the exterior fit is better correlated with the projected string tension. But 
this second choice leads to projected string tensions which are roughly 20% higher 
{(3 = 2.4) than the full string tension. 
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Creutz Ratios at p=2.3 from Products of Z2 Loops 

0.3 I , , , , 




I ' ' ' ' 1 

1 2 3 4 5 6 

I 

Figure 6: Xprod{I,I) at /3 = 2.3. 

A case can be made for either criterion, but the method which we have used 
in the past, now described as a process of "quenched minimization," seems to be 
something of a compromise between the two ahernatives: quenching emphasizes the 
external fit, subsequent minimization the overall fit. The projected string tensions 
obtained in this way also come out roughly in the middle of the two extremes, 
and have good center dominance properties. Apart from center dominance, the 
property of precocious linearity, and the scaling of the vortex density, indicate that 
center projection has isolated the relevant long-range degrees of freedom, rather than 
having them mixed in with physics at all scales. Projected lattices are also found 
to be strongly correlated, in the infrared, among Gribov copies obtained from the 
quenching procedure. 

Nevertheless, the variation in string tension among Gribov copies selected ac- 
cording to different reasonable criteria is much greater than we had expected, and 
the justification for quenched minimization is empirical rather than theoretical. It 
is possible that an improved version of maximal center gauge can be devised which 
retains the appealing "best fit" interpretation, but which softens the contribution 
to the fitting functional at the location of P-vortices. We consider this to be an 
interesting direction for further work. 
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